A note on weakly discontinuous dynamical transitions 



Silvio Franz, 1 Giorgio Parisi, 2 Federico Ricci-Tersenghi, 2 Tommaso Rizzo, 2 and Pierfrancesco Urbani 

1 Laboratoire de Physique Theorique et Modeles Statistiques, 
CNRS et Universite Paris-Sud 11, UMR8626, Bat. 100, 91405 Orsay Cedex, France 
2 Dipartimento di Fisica, Sapienza Universitd di Roma, INFN, 
Sezione di Roma I, IPFC - CNR, P.le A. Moro 2, 1-00185 Roma, Italy 
'Dipartimento di Fisica, Sapienza Universitd di Roma, P.le A. Moro 2, F00185 Roma, Italy 

We analyze Mode Coupling discontinuous transition in the limit of vanishing discontinuity, ap- 
proaching the so called ".A3" point. In these conditions structural relaxation and fluctuations appear 
to have universal form independent from the details of the system. The analysis of this limiting case 
suggests new ways for looking at the Mode Coupling equations in the general case. 



I. INTRODUCTION 



The dynamics of supercooled liquids is characterized by a two step relaxation. After a rapid decay, the dynamical 
correlation function displays a plateau where relaxation is arrested before decaying on a much larger time scale. Mode 
Coupling Theory (MCT) describes the formation of the plateau in terms of a discontinuous dynamical transition 
where the length of the plateau diverges as the temperature become close to the dynamical transition point [l[ . 

The approach and the departure from the plateau are described by power laws, respectively t~ a and t b , where the 
powers a and b are system dependent but obey the universal relation 

r 2 (i + b) _ r 2 (i-a) 

T(l + 2b) T(l -2a)' 1 ' 

The exponent parameter A also appears in replica theory, where it has been related to the ratio between six point 
static correlation functions that can in principle be measured ore computed directly using the Boltzmann measure 
2] . Explicit analytic computations have been performed in mean- field schematic models [JQ and in liquids 0, H| • 

However, the discontinuous glass transitions is not the only possibility. A different transition mechanism is found 
for example in spin glasses with full replica symmetry breaking, where the long time limit of the dynamical correlation 
function passes continuously from zero to a non zero value when the transition is crossed. Within MCT, Bengtzelius, 
Gotze and Sjolander Q have proposed a schematic model whose dynamical transition can be tuned smoothly from a 
discontinuous one to a continuos one through the variation of a parameter. The resulting singularity at the continuity 
point has been named "A3" . 

The study of discontinuous/continuous crossover is not a mere academic exercise. Realistic systems where this is 
found include disordered spin models in presence of a magnetic field, liquids in porous media both in the MCT (ToL ITU 
and in the HNC approximations [H, [l3[ and liquid models with pinned particles [3]. Gotze and Sjogren jl5l| have 
studied the scaling form of the approach to the plateau of the correlation function during this crossover within MCT. 
However, to the best of our knowledge, a full characterization of the correlation function in the a regime below the 
plateau has not been presented in the literature. To fill this gap, we study the properties of MCT equations for weakly 
discontinuous transition. We find that in this case relaxation takes a form which is universal within the theory. In 
the way of this taste we find some new results on schematic MCT equations that have an interest by themselves. 

In Sec. |TT]we set up the problem and discuss the dynamical correlation function for weakly discontinuous transitions 
both in equilibrium and in the aging regime. In section IIIII we extend our analysis to the study of fluctuations and 
compute the four point susceptibility. Finally we draw our conclusions. 



II. MCT EQUATIONS NEAR A CONTINUOUS TRANSITION 



The Mode-Coupling theory postulates that the dynamical correlation and response functions can be obtained solving 
a system of integro-differential equations In the general theory of liquids these are equations for the dynamical 
structure factor and they contain information about the spatial structure of this quantity. However close to the 
transition the spatial structure can be neglected in a first approximation by looking at the peak of the static structure 
factor. This has been first underlined by Bengtzelius et al. but it has been noted also in the framework of the replica 
approach to the glass transition Q . Using this fact one can produce a dynamical equation describing the evolution of 
a single mode, that is called the schematic MCT equation. It is well known that this equation is exactly the one that 
describes the Langevin dynamics of fully connected spherical p-spin model with a IRSB dynamical transition (l6l - [l8| . 
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Approaching the dynamical temperature from above we look at the schematic Mode-Coupling equation for the 
correlation function C(t) 

dC(t) _ ^v^VV^vvi 1 [* a dC '( u 



dt 



= -TC(t) + (1 - C(t))M[C(t)} -^J Q du ^ (M[C{t - u)] - M[C{t)\) , (2) 



where M[C(t)] is the memory kernel that depends on the parameters of the problem 1 , temperature and/or density; the 
initial condition is C(0) = 1. Depending on the nature of the memory kernel, different kind of dynamical transitions 
are possible. 

The simplest scenario, relevant for supercooled liquids is the one of a discontinuous transition. Where above and 
close to the transition the correlation function display a characteristic two step relaxation and ergodicity breaking 
below the transition. 

Above the critical point the asymptotic value of the correlation is given by the unique solution qo of the equation 

q=(l-q)M[q]. (3) 

where we have defined M[q] = M[q]/T 2 . Close to the dynamical temperature Td the correlation function develops a 
long plateau at an intermediate level between 1 and go before relaxing. At the critical temperature the length of the 
plateau diverges. Correspondingly, a second solution q\ > qo to ([3]) discontinuously appears. For this value of the 
correlation one also has 

l = ±(l-q)M[q]\ q=9l (4) 

which expresses a marginal stability condition of dynamics at criticality [l[ . 

Close to the transition, the solution of eq. © in the "alpha regime" , describing the correlation decays below qi , 
verifies the "time-temperature superposition principle" , i.e. it has a scaling form 

C(t,T)wC(i/r (T)) (5) 

where the T a (T) ~ (T — T^) -7 is the relaxation time as a function of the temperature and C(u) is a scaling function 
independent of the temperature. This scaling function can be computed from the equation ([2]) exactly at the critical 
dynamical temperature where one can neglect the time derivative. The MCT equation without the time derivative is 
invariant with respect to rescaling of time t — > at. Accordingly r a is not directly extracted from this approximation 
and one can measure time in arbitrary units (e.g. in units of r Q ). 

Other kind of transitions are possible if the second solution appears continuously. As discussed lengthy by Gotze 
and collaborators [9j, depending on the control parameters in M the transition can change from discontinuous to 
continuous, in passing through a critical point. 

We are interested in the case of weakly discontinuous transitions close to a critical point where q± and go are almost 
degenerate. Thanks to the vicinity to criticality we can characterize these transitions in a universal way. For small 
qi — qo, the exponent parameter A, which is in general determined by the relation 

A = TdM"(qi) 
2(M'( 91 ))3/ 2 



is near to f and both the exponents a and b are close to zero. To the leading order a = b = y ~ (1 — A) ~ y/q[ — qo- 
We choose to parameterize the distance from the critical point by the value of b itself (so that q\ — qa is a vanishing 
function of b in the limit 6—^0). 

As remarked long ago by Gotze and Sjogren, at small argument the function C(i) admits a regular short time series 
expansion in terms of the parameter z — t b , whose coefficient can be computed recursively from Unfortunately 
this expansion is not convergent in the general case, but for b — > we can compute the solution directly from the 
equation. More precisely, we suppose the existence of the limit: 



, lim (C(t,b)-q )/(q 1 -q a ) = G(z), (7) 



1 For example in the p-spin spherical model we have M[q] = pq p 1 /2. 
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with G{z) a well defined function of its argument. 

Let us now rewrite the equation ^ in the a regime. We get: 

C(t) = M\C(t)](l - C(t)) - [ du (M[C(t - u)] - M[C(t)]) . (8) 

Jo du 1 

We now consider the various terms in the equation ([2]). We firstly consider the memory term in the integral; in the 
b — > limit 

M[C(t - u)] - M[C(t)] ~ M'[C(t)](C(t -u)- C{t)) (9) 
cbzM'( qi )( qi -q )^l\n{l-^ 
In an analogous way we have 

C'{u)=^ h -z{ qi -q )^-. (10) 
u dz 

Next we observe that generically at the transition point the function N(C) = —C + (1 — C)M(C) has a single root 
in qo and a double root in qi . For small (ji — go its form should read N(C) — — A(C — <Zo)(<Zi — C) 2 , where by using the 
relations ((3]), (0} and ((U we have A = — ^ 1 1 2^~ A ^ ■ It follows that to the leading order the mode coupling equation 
can be rewritten as 

M'( qi )( qi - qo ) 2 [(l - A)G(1 - Gf - {bzf [G'{z)] 2 [ — ln(l - u)] . (11) 

Jo u 

Now, taking into account that 1 — A = b 2 J Q — ln(l — u)] = b 2 ^- : we obtain the following equation for G: 

G(l - G) 2 = z 1 [G'(z)f . (12) 
This equation is similar to the one found in [l5j and used there to describe the /? regime. Recasting it under the form 

dG dz 



'G{l-G) z 
we find that it admits the solutions 

°«=(^)'- 

The value of zo cannot be computed, as a consequence of scaling invariance of the MCT equation © and we choose 
Zq = 1. We notice that G(z) decreases from 1 to 0, vanishing at finite z — zo- This is not in contradiction with 
the fact that the correlation is positive for all times at finite b, but is a consequence of the fact that we have taken 
the limit b — >• 0. A detailed computation for small but finite b tells us that for z > zq C(z) ~ e ~M z / z o) 1/b _ This 
expression is exponentially small for b — > and corresponds to the simple exponential C(t) ~ e — At/to j n ^ erms f ^ 
where to = z\j h . 

We can compare this asymptotic solution with the Pade approximants of the series expansion of the equation ^fy 

for small values of b. This is done in figure (HU for the schematic F 12 model Q where M(C) = (2A ~ 1 A ) f +C . The 
curves show that the Pade approximants give an accurate description of the function at time smaller that 1, and that 
the limit A — > 1 is achieved smoothly. 



A. Aging 

The previous analysis can be generalized to the aging dynamics. We specialize to the case of the generalized 
spherical p-spin model where the temperature appears explicitly into the equation. The structure of the equation in 
the aging alpha regime is similar to the equilibrium case and one has [19J 

= -TC(t, t') + 0[(ii/'(gi)(l -x)- qo f( qo )x]C(t, (15) 
+Pf'(C(t, t'))(l - qi ) - pf{ qi )(\ - x)C(t, if) 
-px qo f'( qo ) + pxf(C(t, t'))( qi - C(t, t')) 

-frj d* dC Q g ,a) [f'(C(t,s)) - f'(C(t,t'))} . 
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FIG. 1: Scaling function C(z). From top to bottom A = 0.9, 0.99, 1. 
The first two curves are obtained from the (20, 20) Pade approximants of the small time expansion in t b . The last curve is the 

function (yipf j • 



Here /'(C) generalizes the memory kernel M of the equilibrium case. 

The quantity x is the so called fluctuation-dissipation ratio, fixed by the condition that the function 

K(C) = -TC + (3[ qi f'( qi )(l -x)- q f(q )x]C + (16) 
(3f(C)(l - Ql ) - (3f( qi )(l - x)C - pxq f'{q ) - (3xf(C)( qi - C) 

has a double root in C = q±. 

It is well known that equation (|15[) is reparametrization invariant and admit scaling solutions of the form C(t, t') = 
C(g(t) — g(t')) where the reparametrization function g(t) is left undetermined. The short time expansion of the 
equation predicts a behavior of the kind 



C(u) = q x + (u) b , (17) 



where b is determined by the condition [l£ 



T f'"( qi ) Tjl + bf 

2//'( ?1 )i ~ r(l+26)' U8j 

As in the equilibrium case, for q\ close to (jo the function K(C) behaves as K{C) — A(C — qo)(qi ~ C) 2 ). We 
can suppose that C becomes an analytic function of z = (g(t) — g(t')) b . Notice that if the function g(t) is such that 

g" {t) / g' (t) 2 « 1 for large t, then one can equivalently write z = ■ ^ e can then define the scaling function 

G(z)= hm ^'^-go (1Q) 
b-Ki, t.f-foo qi-q 

z = (g(t)-g(t')) b 

and repeat verbatim the analysis of the equilibrium case. It turns out that the equation verified by G coincide with 
the one found at the critical point. A fortiori, the same is true for the function G(z). 



III. FLUCTUATIONS 



In this section we would like exploit our analysis to investigate fluctuations in the alpha regime. In the last years, 
research has concentrated in the study of fluctuations of the time dependent correlation functions in terms of 4-point 
functions. As often in disordered systems one can define different kinds of correlation functions with a-priori different 
scaling properties. It has been recently proposed that it is useful to disentangle the fluctuations of correlations with 
respect to thermal noise for fixed initial condition from the fluctuations with respect to initial conditions [20j | . 
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Denoting by (•) the thermal average for fixed initial condition (iso-configurational average) and by [•] the average 
initial condition, we define [2(| [2l[ 

Xt h {t) = [{C{tY)]-[{C{t))*], (20) 
Xh et(t) = [(C(t)) 2 }-[(C(t))} 2 - 

A theory for this kind of fluctuations in the beta regime has been proposed in [2p| . using a " reparametrization 
invariant" formulation where time is eliminated in favour of the average correlation function. Within a gaussian 
fluctuation theory it is found that the singularity of Xhet doubles the one of Xth- 

The basic observation allowing to study now the functions in the a regime is the fact that, as proposed in (2^, the 
leading behavior of Xth{t) can be obtained as 

/ \ dC(t) 

Xth{t) oc -^.Z. (21) 

Before exploiting this relation we would like to note that it appears naturally in the theory put forward in [2(J. In 
that context that fluctuations can be described through a field theory where the correlation function, which plays the 
role of fundamental field, couples linearly to the temperature. Moreover, the dependence with respect to the initial 
configuration turns out to be parameterized by a random variation of the temperature. This has the consequence 
that the susceptibility Xhet is the square of the thermal one multiplied by the variance of the random temperature. 

While these consideration strictly hold for the beta and early alpha regime, the time-temperature superposition 
principle shows how the correlation is very sensitive to any temperature change which can induce large changes in the 
relaxation time. This is a sort of "beta imprinting" indicating that large fluctuations of the correlation fluctuations in 
the alpha regime could be just consequence of fluctuations in the initial time of relaxation. In last instance this is a 
consequence of the emerging scale invariance of the MCT equation when the critical temperature is approached. We 
see here a link with the theory of fluctuations during aging dynamics below T d developed by Cugliandolo, Chamon 
and collaborators [23-26] where fluctuations are ascribed to the large time emergence of reparametrization invariance 

With all this in mind, we can write: 

Nxtk{t) . fS^Q) , (22) 

N X het(t) = [5T 2 ]xth{t) 2 ■ 
Using the relation r(T) ~ (T - T d ) with a « b for b -> and C(u) = G(u b ), one gets 

xt h (t) = 2 Y^Y d iqi - *^G'(z) - 2 Y^Y d iqi ~ qo) ^ {l ~ G) ' (23) 

Xhet{t) = 4[<5T 2 ] {T \ d)2 (<li - Qo) 2 G(l - G) 2 . 

The divergence as a function of T — T d which just depends on the power law behavior of the relaxation time, confirms 
the direct dynamical analysis of 22]. 



Notice that for a finite system the divergence should be cut-off by a function of the volume. It was found in [27 1 
and (2(| that the scaling variable describing the cross-over is x = (T — T^N 1 / 2 . This predict an alpha relaxation 
scaling at T d where Xth ~ -j= and a finite Xhet- 

In 20] it was shown that if C(t) follows a bimodal distribution as it would be implied by a simple jump process, 
one should expect the dependence Xhet ~ G(l — G). Notice the form we find differ form this expectation. 

We would like to remark that while the square root behavior of Xth at small G is only valid in the limit of small 
Si ~ Qo that we are considering, the linear behavior for G ~ 1 is more general: it is a consequence of the initial power 
law relaxation of the correlation function C(t) = q± — at , that hold whenever there is a discontinuous transition. As far 
as the small C behavior for finite b is concerned, the final exponential relaxation suggests a behavior xth ~ — ClogC. 



IV. CONCLUSIONS 



The point where the discontinuous transition becomes continuous can be seen as a critical point for Mode Coupling 
Theory. As such universal properties emerge which do not depend of the details of the model [l|. In this note we 
have computed the scaling functions for the correlation function both at the MCT transition and in the aging regime, 
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finding that they take the same universal form. We have also analyzed the behavior of fluctuations, finding general 
expressions of the four point functions as a function of the correlations. 
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